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Abstract 

We consider a system of diffusing particles on the real line in a quadratic 
external potential and with repulsive electrostatic interaction. The em- 
pirical measure process is known to converge weakly to a deterministic 
measure-valued process as the number of particles tends to infinity. Pro- 
vided the initial fluctuations are small, the rescaled linear statistics of the 
empirical measure process converge in distribution to a Gaussian limit 
for sufficiently smooth test functions. We derive explicit general formulae 
for the mean and covariance in this central limit theorem by analyzing a 
partial differential equation characterizing the limiting fluctuations. 



1 Introduction 

We consider the following system of n ltd equations: 

d\l = -^=dBl - Xjdt + — V — ^— , for i = 1, . . . ,n. (1) 

Here are independent, standard Brownian motions and a and (3 > are 

real parameters. These equations model the dynamics of n diffusing particles 
on the real line with a logarithmic interaction potential, u{x) = — h log \x\, con- 
strained by a quadratic external potential v n (x) — at inverse temperature 
P. Cepa and Lepingle proved that the order of the particles is almost surely 
preserved for all times t > 0. The stationary solution to (QJ has distribution 

— exp i -j3 I ^ v n (Xj) + W ( A * ~ A j) f II dXl 

=i n i^-^i^xp(-^E^}n^' ( 2 ) 

At l<i<j<n y j=l J i=l 

where is a normalizing constant (the partition function) and dX denotes 
Lebesgue measure. 
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For the specific parameter values = 1,2 and 4 this model can also be 
interpreted in terms of matrix valued stochastic processes (Dyson Brownian 
motion). Let A4 n (0) be the set of all n x n real (0 = 1), complex (0 = 2) and 
quaternion (0 — 4) matrices respectively and S n (0) the set of self dual (with 
respect to conjugate transposition) elements in M n (0)- The Gaussian Orthog- 
onal (0 = 1), Unitary (0 — 2) and Symplectic (0 = 4) ensembles, GXpE n (a 2 ), 
with Xp = 0,U,S for = 1, 2, 4, are the probability distributions 

d^ n (M) = — exp |-|^TrM 2 ] dM 

on S n (0), where dM = n"=i * ni<i<,-< B ' ' ' dM iP is P roduct Lebesgue 

measure on the essentially different members of M = (M-j, . . . , M^- Let 
M t = {Ml 1 ) ij be an S„(/3) -valued Ornstein-Uhlenbeck process, i.e. satisfying 
the SDE 

dM t = -M t dt + -%==d(B t + B*), 
V0n 

where B ( isannxn matrix, the elements of which are independent standard real 
(0 = 1), complex (0 = 2) or quaternion (0 = 4) Brownian motions and i? t * is the 
conjugate transpose of B t . Then the eigenvalues {Aj}™ =1 of M t satisfy (see 
|E]). For instance, if M £ GXpE n (al) we will have M t £ GX(3E n {e~ 2t {al - 
a 2 ) + a 2 ) for all £ > and if M £ S n (0) is fixed, equation © has solution 
M t = e-'Mo + N t , where N t £ GX fj E n (a 2 (l - e" 24 )). 
We define the empirical measure process 

i=l 

To capture the asymptotic properties of the model on a global scale as n — > 
oo, one is interested in studying the limiting behaviour of the linear statistics 
(X^ 1 , f) = \ Y!,?=i f(K) where / is a bounded continuous real test function. 

Example Define a deformed GUE to be an ensemble of Hermitian matrices 
M a = M + D a where M is distributed according to the GUE% n {X) and D a — 
(dy)f"_j is a fixed 2n x 2n diagonal matrix with 

. (a for 1 < i < n 
11 \ —a for n + 1 < i < In 

Then the eigenvalues of the rescaled matrix M a /\/a 2 + a~ 2 correspond to the 
particles in our model with initial distribution Xq = X^ n = + Si) at 

time t = log yj\ + (acr)~ 2 . The local behaviour of the eigenvalues in this model 
have been studied in and [2] and it is known that the limiting eigenvalue 
density of M a as n — > oo is supported on two disjoint intervals if < a < 1 
and on one single interval if a > 1. In other words, suppX t grows from the two 
starting points {—1, 1} at time t = into two disjoint intervals that join at time 
t = lo gv /l + (a<7)- 2 . 
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In the stationary case, it is a classical result that (X n , /) converges in dis- 
tribution to L /d/i, where 



is the Wigner semi-circle law. More generally, for any initial asymptotic dis- 
tribution of particles, Xq, the asymptotic particle distribution X t at each time 
t > is uniquely determined by Xq and converges weakly to [i as t — > oo (see 
Theorem 12 .11 for a more precise statement). 

A natural question is if there is a limiting distribution of the rescaled linear 
statistics. For ease of notation, we introduce the fluctuation process Y™ = 
n(X"—X t ), which takes signed Borel measures on R as values. We are interested 
in the limiting distribution of the random variables 



where / is a test function from an appropriate class, as n — ► oo. Note that there 
is no normalization of the linear statistics here; this reflects the very regular 
spacing of the particles and is typical of related models. 

Israelsson Jl] showed that ((F t ™, /i), . . . , (Yj™, /&)) has a Gaussian limit, 
provided the the test functions fj are 6 times continuously differentiable and the 
initial distributions Xq converge sufficiently fast to Xq, Although establishing 
existence and uniqueness, he does not characterize the limiting distribution- 
valued Gaussian process Y t very explicitly. In this work we derive explicit general 
formulae for the mean and covariance of the finite dimensional distributions of Y t 
by analyzing the partial differential equation arising in Israelsson's proof. These 
formulae generalize many similar results obtained for various special cases of 
our model by completely different methods, some of these are briefly discussed 
below. In particular it is worth noting that our results hold for all values of 
the inverse temperature (3 and in the non-equilibrium case with arbitrary initial 
particle distribution Xq. 

Most of the previous related results pertain to specific matrix models and are 
restricted to the cases f3 = 1 or f3 = 2, The asymptotic global fluctuations for 
various ensembles of Hermitian and real symmetric matrices have been exten- 
sively studied, see e.g. [Hj, [E], [Tfl| - [15] and JZj- In a recent paper [Tjj. Bai 
and Yao consider N x N matrices with zero mean, independent, not necessarily 
identically distributed entries such that the diagonal elements all have the same 
variance a 2 /N and the off-diagonal elements have variance 1/N (real symmetric 
case) or uncorrelated real and imaginary parts each of variance 1/2N (Hermi- 
tian case). Such models are known as Wigner ensembles. Under some fourth 
moment conditions, they provide a central limit theorem and give explicit mean 
and covariance formulae, which agree with those of Corollary 12.411 with (3 = 2 
and At = 0. Under the assumption of finite moments of all orders of all matrix 
elements, a more general class of ensembles of symmetric matrices is considered 
by Anderson and Zeitouni in [lj. Here the variances of all matrix elements, and 
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the means of the diagonal entries, are allowed to depend on position. Spohn 
|18j derives an expression for the covariance of the Gaussian fluctuations of our 
model in the hydrodynamic limit, but deals only with the case (3 = 2 and (time 
dependent) equilibrium fluctuations. 

The few previous results available on the general (3 case are restricted to an 
equilibrium situation. For the corresponding model on the circle, Spohn 
can handle the general (3 case, again in the hydrodynamic limit at equilibrium. 
By expressing the equilibrium model in terms of ensembles of tridiagonal real 
matrices, Dumitriu and Edelman [2] are able to find the general (3 global fluctu- 
ations for polynomial test functions (this corresponds to the leading order term 
in Proposition 123- Johansson \12\ considers a more general model correspond- 
ing to the equilibrium measure <|2J but with the quadratic external potential 
v n replaced by a general polynomial of even degree and with positive leading 
coefficient. For the case of quadratic v n his mean and covariance formulae agree 
with the fixed t equilibrium case of the model we discuss. In Johansson's model 
the variance is universal in the sense that it does not depend on the details of 
the potential, provided the support of the equilibrium measure is a single inter- 
val. In our model however, the variance at every finite t depends on the initial 
conditions (see Pronosition l2.6l) . For instance, even though the eigenvalue den- 
sity in the deformed GUE example will in finite time be supported on a single 
interval, the fluctuations remember the initial particle distribution for all t > 0. 
Thus the time evolution of the variance is determined by the geometry of the 
initial distribution; this structure is reminiscent of the role played by the bound- 
ary conditions in determining the fluctuations of the height function in discrete 
plane tiling models such as Kenyon's There the fluctuations converge to a 
Gaussian free field for a conformal structure determined by the boundary. 



Acknowledgment I wish to thank Kurt Johansson for invaluable advice 
throughout the preparation of this work. 



2 Main results 

In order to formulate our results we need the following Theorem, referred to in 
the introduction: 

Theorem 2.1 (Rogers and Shi, Cepa and Lepingle, 0) Suppose that Xq 
converges weakly in M., the space of Borel probability measures on R with the 
weak topology, to a point mass Xo at an arbitrary element of M. . Then there is 
a family {Xt}t>o C M., depending only on Xo and converging weakly as t — ► oo 
to the Wigner semi-circle law, [i, such that for each t > 0, X™ converges weakly 
to X t in M. as n — > oo. X t is uniquely characterized by the property that its 
Stieltjes transform, 

M = M(t, z) = [ dXt{x \ (t, z) € [0, oo) x (C \ R), 



4 



solves the initial value problem 

M t = {2a 2 M + z)M z + M, t>0 

M (0,z) = f^^M. (5) 

We fix some terminology that will be used throughout the rest of this paper. 
Let X be a given Borel probability measure on M and define f2 = C \ K. 
Put f(w) — J ^x-w > w ^' f wn ^ ^ e a holomorphic function. It follows 
from Theorem 12.11 that for every t > 0, M(t, •) = J ^j[£lg) j g a well-defined 
holomorphic function in Q, so we can define a family {/it}t>o of holomorphic 
maps in Q, 

h t (z) = ze f + er 2 (e* - e _ *)M(t, z). (6) 

Proposition 2.2 For euen/ t > 0, /it(f2) C f2 and i/ie relation gt ° h t = id 
holds, where 

g t (w)=e~ t w-a 2 (e t -e- t )f(w). (7) 
Define = g t2 ° /i tl /or ti > t 2 > 0. T/ien /i^ = /it 2 ° ^ti ■ 

Proof This is a step in the proof of Theorem 12.31 | 

Recall the definition of the Schwarzian derivative: 

Definition Let v be a univalent function in some domain of the complex plane. 
The Schwarzian derivative Sv of v is defined as 



(Sv)(z) 



v"'(z) 3/t/'(z) x2 



v'(z) 2 \v'{z) 



We define the generalized Schwarzian derivative, also denoted Sv, to be the 
following function of two complex variables: 



{Sv){z 1 ,z 2 ) 



dz 1 dz 2 ° I 21—^2 / (tl(zi)-D(z 2 )) 2 (Zl— Z 2 ) 

(8) 



Definition A probability measure on the space 5' of tempered distributions 
is called an initial measure if its law is a point mass at a v 6 5' such that 

(M)=o. 

We can now state the main result, giving expressions for the mean and 
covariance of the finite dimensional distributions of the Stieltjes transform of 
the limiting fluctuation process Y t . 
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Theorem 2.3 (Mean and covariance formulae) LetY™ = n(X™—X t ) where X" 
is the empirical measure process and X t its weak limit. Suppose Yq converges 
weakly in S 1 to an initial measure Yq and that there is a constant C such that 
for every n and z = a + bi, b ^ 0, the it, 

E 



dY n (x) 2 



x — z 



C 



holds. Then Y" converges weakly to a Gaussian distribution-valued process Y t 
(see Theorem, I3.il for the full statement), and for < tk < ife-i • ■ ■ < t\ and 
z = (zi, . . . , Zfe) € (C \ R) h the Gaussian random vector 

U=(U 1 ,..., U k ), where Uj = (Y tj , -J— ) 



has mean 



1/2 \ h'lizj) ti t .( Zj ) 



(9) 



and covariance matrix 

a a 2 a 2 ( h tj { Zj ) - h tl ( Zl y 

2 



In particular, 



Proof See sectional 



p h% [z s ){Sht t ) ifl>j. (10) 



Var(C^) = -L (Sh t .) (zj). (11) 



Remark An interesting consequence of equation ltTT|l is that a complete knowl- 
edge of the variance of the Stieltjes transform of Y t as a function of z in the 
upper (or lower) half-plane at any fixed time t > 0, will uniquely determine the 
initial particle distribution Xg: By (1 1 H this function is the Schwarzian deriva- 
tive of some analytic function, which is unique up to composition with an arbi- 
trary Mobius transformation (this is a well-known property of the Schwarzian 
derivative). It is easy to see that this determines h t uniquely, and this in turn 
determines the Stieltjes transform / of the initial particle distribution X . 

Remark Another model, concerning eigenvalues of non-Hermitian complex ma- 
trices, where a similar variance formula involving the Schwarzian derivative oc- 
curs is studied in 20J. 

Let 

//*(*) = 
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denote the Stieltjes transform of the Wigner semi-circle law, fi. (Here \f means 
the branch of the square root for which ^(^fz) > iff 3z > 0, defined for 

z e C\ (-oo,0).) 

Corollary 2.4 (Equilibrium fluctuations) Let z 1} z 2 € f2 and At > be given. 
Put ti = t + At and t 2 = t. Under the same hypotheses as in 

Theorem\EM the 

asymptotic mean and covariance as t — ► oo are given by 
and 



lim A12 

t — >oo 




f3(^U( Zl )U(z 2 )e-^~l)^ { ->> 



Proof This is just a calculation using Theorem 12.31 and the fact that X t con- 
verges weakly to the semi-circle law fTheorem l2,lll . § 

The previous results can be expressed in terms of integral formulae for the 
fluctuation process acting on analytic test functions. Suppose that for each t > 
there is a compact set Ct C K such that suppX t C Ct and, with probability 1, 
supply C Ct; we say that the fluctuation process Y t is compactly supported. Let 
t\ > t 2 > and for i = 1,2 let ji be a closed simple curve in the complex plane, 
the interior of which contains C^ , and let Di be a simply connected domain 
containing 7j, 

Theorem 2.5 (Integral representation) Suppose thatY t is compactly supported. 
Let F\ and F 2 be analytic in the domains D\ and D 2 defined above, respectively. 
Define the random variables Z\ = (Yt t ,Fi) and Z 2 = (Yt 2 ,F 2 ). Then 

Cov(Z l7 Z 2 ) = <j) j> {F 1 (g tl (w 1 )) - F 2 (g t2 (w 2 ))) 2 (Sg t2 ) (w 1 ,w 2 )dw 2 dw 1 , 

(14) 

where Ti — /it ( (7i)- For Z\ = Z 2 this reduces further to 
VarlW,,^), ' / f ^.te,K))-F,b„W) y^ droi . (15) 

Remark Note that the covariance depends only on the initial distribution Xq. 
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Proof Representing F\ and F 2 by the Cauchy integral formula as contour in- 
tegrals along 71 and 72 we can use the linearity of Y t and Fubini's theorem to 
obtain 



Cov(Zi,Z 2 ) = E(ZiZ 2 ) - EZ1EZ2 



E 



y J_ / F 1 (z)dz j_ r 



F 2 (z)dz. 



x — z 



-E 



1 



F 1 {z 1 )dz. 



27TJ ./ 7l 



E 



(Yt 2 , 



F 2 {z)dz^ 



2ni ./ 72 f - : 



-E 



F 1 (z 1 )F 2 (z 2 ) 

^z7~- 

Fi(zi)F 2 (z2) 



E 



E 



1 



(27Ti) 



«2 ' 

z 2 - 
k\ 2 dz 2 dz\, 



dz 2 dz\ 



(16) 



where 



A = 2 ^ 2 1q 1 ~ h t2 (z 2 ) 

(3 dzidz 2 



Kl(z 1 )-z 2 J 

by Theorem l2.3l Since, for fixed z 2 , A% 2 is the derivative of an analytic function 
of z\ in a domain containing 71 , we note that 



{F 2 {z 2 )f k 12 dz 2 dz x 



71 J 72 



1 



(F 2 (z 2 )) 2 </, _ _ 

72 J^i l^tx^i) - h t2 (z 2 ) h\\{z x )-z 2 



dz\dz 2 = 0. 



Similarly, 



<f <f (F 1 (z 1 )) 2 A 12 dz 1 dz 2 = 0, 



so we may substitute — \ (F\{z\) — F 2 (z 2 )) 2 for the factor -^1(21)^2(22) in equa- 
tion iflfil) . which gives 



Cov(Zi,Z 2 ) = 



1 



2 I ti tl (zi)K 2 (z2) 



-F 2 (z 2 )Y 



{h tl (zi) - h t2 {z 2 )f {z 2 -h t t \{z 1 )f 



dz 2 dz\. 



In the variance case, Z\ = Z 2 , this reduces further since the second term of the 
integral becomes 



Fi(2i)-Fi(2 2 ) 

z 2 - z\ 



dz 2 dz\, 
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which vanishes by the analyticity of F±. With the change of variables w\ = 
h tl (zi), w 2 = h t2 {z2) we arrive at the expression l|14j) , | 

As an application of Theorem 12. 51 we can consider polynomial test functions. 



Proposition 2.6 (Variance for polynomial test functions) Suppose Y t is com- 
pactly supported. Then for n — 1,2,... 



4a'" \r^n/z 

Var((Y t ,x n )) = < ^s=x°\ n /2+s) 



e t R n (t,Xo) if n is even 



(17) 



if n is odd, 

where Ru is bounded in t and depends only on moments of Xq up to order 2k— 1. 

Proof Let F%(x) = i^Oc) = % n and t\ = t% = t. After the change of variables 
Zi = e* /wi, formula l(T5l) in this case reads 



Var((Y t ,x n )) 



4tt 2 /3 J M=r 



|*l|=r J\z 2 \=r 



g t (e t /z 1 ) n -g t {e t /z 2 y 
Zl - z 2 



dz\dz2, (18) 



where r > is such that suppAo C (—e t /r,e t /r). By definition of gt, we can 
expand gt(e* / ' z) in a Laurent series, 



(«7*) = E 



a k z 



fc=-i 



where o_i = 1, a = and a fc = <7 2 (1 - e -2«) e -(fc-i)t J x k-idX (x) for fc > 1, 
so equation i|18H can be written 



Var((y t ,x"» 
1 



r J \z 2 \=r 



2_j a fei • • • a fc„ 



^felH hfen ,felH hfcn 



- ^2 



dz\dz2- 
(19) 



For given integers K and J a simple combinatorial argument and the residue 
Theorem show that 



1 



'X y -2 



dz\dz2 



I \ Zl \=r J \z 2 \=r \ z \- z 2 J \Z\- Z 2 

This means that we can write equation i|19ll in the form 



\K\ HK=-J 
otherwise. 



Var((Y t) x n )) = - V M_ s>n A s , n , 



(20) 



9 



where 



Asm — flfei • • • a-k„ , 



&! + ■■• + k n = S 

h > -l 

is a finite sum with the following structure: Since in the limit t — ► oo, a\ — > er 2 
and afe — > for fc > 1, all terms contributing to v4 Si „ tend to exponentially in 
t unless ki = ±1 for i = 1, . . . , n. If n + s is odd there are no such terms, and if 
n + s is even there are (n+s) choices of (ki, . . . , k n ). Thus 

2 ' 

( if n + s is odd 

lim A s n = S I n \ n+s -r i 

t^oo s ' n [n+s) a ^ if n + s is even. 
Inserting this into equation l|20"jl gives ifTTjl . | 

Remark The t — > oo limit in this formula agrees with the variance formula of 
Dumitriu and Edelman [2] and that of Johansson ,12] which asserts that 

1 00 /2 \ 2 
Var((Yt, /i(a;))) = — V - / /i(2<7cos(0)) cos(fc6)de 
2/3^ Wo / 

for an appropriate class of real test functions h. Indeed, rewriting h(x) = x 2n 
in terms of Chebyshev polynomials, the asymptotic variance is recovered from 
Johansson's result. 



3 Proof of Theorem 2.3 

The proof of Theorem 12.31 relies on the characterization of Y t provided in Jl] 
to prove existence and uniqueness of this process. For convenient reference we 
restate this result. 

Theorem 3.1 (Israelsson, jH]) Suppose that the sequence of measure-valued 
random variables Yq converges weakly in S 1 to an initial measure Yq. Suppose 
further that there is a constant C such that for every n and z = a + hi, b ^ 0, 
the inequality 




holds. 

Then Y t n converges weakly to a Gaussian distribution-valued process Y t in the 
sense that for any 6 times continuously differentiable, rapidly decreasing real test 
functions fj, the random vector (J fi(x)dY£(x), . . . ,J fk(x)dY t ™(x)) converges 
in distribution to ( (Y tl , f %),... , (Y tk , /&) ) . Furthermore, the convergence extends 
to test functions of the form z € O = C \ K, and Y t is uniquely character- 
ized by its action on such functions by the following property: Let < t m +k < 
tm+k-i ■ ■ ■ < Un+i < h < T be given and for s = (si, . . . , s m , . . . , s m +k) G 
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C m+fe , z = (zi, . . . , Z m , . . . , z m+ k) € £l m+fe and £ m +i < £ < £1 define the func- 
tion 



(f)(t 7 Si 7 . . . , S m , Zi , . . . , Z m ) 



E 



I 3=171+1 

Then <fi satisfies the PDE 



3 = 1 



Of 



E 

3=1 



g2 aM(M3) \|j 



+ ( Zj +2a 2 M{t, Zj )) 



d<f> 
dzj 



SjdX t (x) 2a 



SjSidXtjx) 

t (x-zrf (3 j~[^-[J {x - zj) 2 (x - z t f 



(21) 



Remark This is a slight reformulation of Israelsson's result: He allows for Yq 
to be random and works with real and imaginary parts of the complex functions 
in order to ensure that the characteristic function be a priori well defined. 
However, once it is established that the distributions are Gaussian for such test 
functions if Yq is an initial measure, 4> will be a well defined entire function of 
s for test functions ze!l. The argument leading to equation il2l|l is then 
identical to that in Israelsson's proof, but this form is convenient for finding 
explicit solutions. 

Remark There is a numerical mistake in Israelsson's derivation of equation l(2Tl) 
which has been corrected here; all occurrences of the factor § in the equations 
on page 51 and onward in pi] should be replaced by a. 

Israelsson's method is similar to, although technically more involved than, 
that used by Rogers and Shi to prove Theorem 12.11 The proof deals first with 
convergence for test functions in S and the convergence is then extended to the 
larger class of test functions. The existence of a limiting process Y t is shown by 
first establishing that the set of probability measures associated with the family 
P"t n }nLi 1S tight in C([0, T],S'). Then it is shown that any subsequential limit 
Y t must satisfy a certain martingale problem. By tightness, weak convergence of 
all finite dimensional distributions will entail the existence of a unique weak limit 
Y t . This convergence is established by proving that any subsequential weak limit 
of the sequence Y" has the property that for every / € S and t € [0, T] there 
is an approximating sequence {g^^i of linear combinations of functions of the 
form 73^ such that (Y t ,gj) converges weakly to (Y t ,f). This gives rise to the 

PDE ll2lll for the characteristic function of the Stieltjes transform of Y u which 
is shown to have a unique solution. Given an initial point mass distribution for 
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Yo, the resulting initial value problem is shown to have a unique solution which 
is the characteristic function of a Gaussian vector. Hence the finite dimensional 
distributions converge weakly by the approximation property and the theorem 
is proven. 

By solving equation (|2ll) under appropriate initial conditions, we will be 
able to find expressions for the mean and covariance of the finite dimensional 
distributions of Y t . 

Lemma 3.2 For any fixed to > 0, let (f>t (s,z) be a given analytic function 
defined for s = (s±, . . . , s&) € C k and z = [zi, . . . , Zk) £ Q k and let U = 
{(t,s,z)\t > t ,s € C k ,ze n k } and V = {(t ,s,z)\s G C fc ,z S fl k } C dU. The 
initial value problem 



dt 

3 = 1 



g[„(i + ^2!^fii)^ + fe + ^( ( , J9 )) 

n- 2^ 2 -A [ SjdX t (x) 2cr 2 J^J^ f SjSidX t {x) 

21(7 (^F"1T§§ 



(x - Zj) 2 (x - z{f 



4> in U , 



<f>(t ,s,z) = 4>t (s,z) onT, (22) 
has the following unique solution: 



k k 



(f)(t,s,z) = (f) t0 ^■h t t °'(z),h t t °(z)^exph^2 Sj ^- 2J2J2 s J s l A 3l > , (23) 



3 = 1 3 = 1 1=1 



I (2 ^ 



where 

Aji = A,j = |(,S/i?)fe,zi), (25) 

and s ■ hi (z), hl°(z) is shorthand notation for (sihl (zi), . . . Sfc/i'° (zk)) and 
(h t °(zx), . . . ,h°(zk)) respectively. 

Proof The equation is linear and can be solved with the method of charac- 
teristics. Fix (t, s,z) € U and let 4>{t) = (f>(x(r)) be the solution along the 
characteristic x(t) — (t(r), s(t), z(t)) through that point. By J22J, the equa- 
tions for x(t), if we choose t(r) — r, become 

-^2 = - Zj{T )-2a 2 M(r,z j (r)) (26) 

and 



dr V dzj 



= - Sj (r) 1 + 2^ — (T, Zj (r)) , (27) 
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while the solution 4>{t) = 4>(x{t)) along the characteristic is given by the equa- 
tion 

#(r) f 3 (2 \A f Sj (r)dX t (x) 



1 V/3 z^jy (z-^-(t)) 3 



It will be convenient to solve these equations for all r > and impose the 
initial condition at the end. It may seem difficult to find solutions in closed 
form because of the dependence on the evolution of Xt, which is known only 
through the property of having Stieltjes transform satisfying (JjJ- As we will 
show however, all dependence on X t can be expressed in terms of M and, more 
crucially, the evolution of M along the characteristic is particularly simple. For 
the first point, it is easy to see by algebraic manipulations that 

= ^M zz (t, Zj ), 



dX t (x) 


i d 2 f 


r dx t ( x ) 


(X - Zj) 3 


~ 2 dz] { 


J {X-Zj) 


dX t (x) 


1 d 3 f 


r dx t { x ) ' 




~ 6 dz 3 t 





M zzz (t, Zj ), 



and, with a little more effort, 

dX t (x) _ {2(M(t, z 3 ) - M(t, zi)) M z (t, Zj) + M z (t, z t ) 



(x- zj) 2 (x - zi) 2 \ {zj-zif (Zj - zi) 2 

if Zj 7^ zi. (Differentiating under the integral sign is clearly justified here since 
all integrands are bounded.) Assuming without loss of generality that Zj ^ zi 
if j ^ I, equation l|28j) can thus be written 

1 #(r) 



(r) dr 



2 - 2 ( \ - E ^M^tMt)) S -^M zzz (r,z 3 (r)) 

3=1 \3=1 

2(M (r, Zj (r)) - M(r, zj(t))) M,(t, zj(t)) + M z (t, z,(t)) 



fe(r)-^(r)) 3 (*i(r)-*i(r)) s 



(29) 



The equations Il2(it and 112 711 can now be integrated with the aid of © defining 
the evolution of M (r, Zj(t)). Fix Zj = z and put M(r) ee M (t, z(t)), M z (t) ee 
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dM (g'*( T ^ and so on for all partial derivatives of M{t,z). Differentiating, we 
have by the chain rule and equations and Il2(it : 

= M z {T)z'{ T )+M t {T) = M z (t){z'{t)+2<j 2 M(t)+z{t))+M{t) = M(r). 

CIT 

or in integrated form simply 

M(t) = M{t,z)e T - t . (30) 
With I.'IOII substituted into 12611 . the latter equation can be integrated to yield 

z(0) = ze* + cr 2 (e* - e"*)M(t, z) = h t (z). 
Using this initial condition, equations (l26l) and l(3"0j) give the explicit expression 

z(t) = e- T z(0) - a 2 (e T - e^)/(z(0)) = g T (z(0)) = g T (h t (z)) (31) 

for the characteristic. In particular, taking r = t gives z — g t (h t (z)), and 
since there is a unique characteristic through each point in U it follows that for 
t > t\ > to, h t (z) — h tl (gt 1 (ht(z))), which is the assertion of Proposition 12.21 
Note that this provides a method of calculating h t (and M(t, z)) by finding an 
inverse of the explicitly defined function g t . 

Since the function (t,z) >— > gt{z) — er l z — <r 2 (e* — e^ t )f(z) is C°°, it fol- 
lows from implicit differentiaton of the relation g t (ht(z)) — z that the order of 
differentiation can be interchanged in the mixed partial derivatives of h t (z), in 

particular 9 d ^ z k Z ^ = ° g^dt^ ^ or ^ ~ 1, 2, 3. Using this and differentating 
equation © gives: 

M zt = M tz = (2a 2 M z + l)M z + (2a 2 M + z)M zz + M z , 
M zzt = M tzz = (6a 2 M z + 3)M ZZ + (2a 2 M + z)M zzz , 

and 

M zzzt = M tzzz = 6<j 2 {M zz ) 2 + {8cr 2 M z + A)M ZZZ + (2c 2 M + z)M zzzz . 

With these equations we can obtain ODE's for M z , M zz and M zzz in a com- 
pletely analogous fashion: 

= M zz {t)z'{t) + M zt (r) = 2(a 2 M z (r) + l)M z (r), (32) 
= M zzz {t)z'{t) + M zzt {r) = (6ct 2 M z (t) + 3)M„(r), (33) 



dr 
and 

dMzzz(T) = M zzzz {r)z'{T)+M zzzt {r) = &a 2 {M zz {T)) 2 +A{2a 2 M z {T)+l)M zzz {r). 
dr 

(34) 
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Putting w = h t (z), equation 1-iOII can be expressed 



M(r) = f(w)e\ 
and equations (l32l) through l(34*j) can be integrated to produce 



M z [t) =f'(w) - 
M zz (t) 



9t{ W ) 

f"(w) 



(g'rM) 3 



and 



M zzz (t) = 



/'» + 3(/"H) 2 ^ (r ' " ' ] 



9r( w ) 



(35) 

(36) 
(37) 
(38) 



where g T (w) = e T w — a 2 (e T — e T )f(w). Inserting into equation (1271) and 
integrating we get 

We can now finally express the right hand side of equation Q29I1 as an ex- 
plicit function of r by plugging in our expressions ll3"T)l and lj35l) through l(39l) 
derived for the evolution of sj, Zj and z-derivatives of M along the characteristic. 
Integrating we see that 



where 



and 
III 



2a 



log 

I = 2ia 2 
II 

2 ft 
I 

to 



</>(t,B,z) 



4>{t ,s(t ),z(t )) 

t k 

I 

o j= i 



I + 11 + III, 



2a_ 



l )j t Y. Sj ^M zz {T, Zj {T))dT, 



to 



E 



6 



-M zzz (t, Zj{r)) dr, 



(40) 

(41) 
(42) 



2(M(T,Zj(T))-M(T, Zl (T))) 

M z (t, z 3 -(r)) + M z (t, z,(t)) 
{z 3 {T)~ Zl {T)f 



dr. (43) 



To calculate these integrals we first note some immediate consequences of the 
definitions of the function h^ = g to o h t and the generalized Schwarzian deriva- 
tive: 

hl°"(z)_ 1 fg'lW g»{w)\ , 44) 



hl a '{z) g'tH liW sIH/' 
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and 



(Shf)(zi,Z2) 



9t{wi)g' t {w 2 ) 



{{Sgt Q ){w 1 ,w 2 ) - {Sg t ){w u w 2 )) ■ 



(45) 



where Wi — h t (zi). Using the change of variables x = 6(r) = cr 2 (e 2r — 1) we can 
now calculate the integrals on the right hand side of equation J32J). First, we 
note that by lEPfjl . 



s{t)M zz {t)cLt = s 
sf"(w) [W 



t „/ 



dr 



t0 g' t (w) (g' T (w)r 
dx sf"H 



2a 2 g' t (w) J b{to) (1 - xf'{w)Y 2^g' t (w)f(w) 



1 - xf'(w) 



b(t) 



fc(*o) 



e-%t)f"(w) 



s (z) 



2^g' t {w)\e- t (l-b{t)f'{w)) e -*o(l - b(t )f'(w))J 2a 2 {z) 
This means that 



2ia 2 







0E- 



h\°"{z 3 ) 



(46) 
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For the second integral, 



S {T) 2 M zzz ( T )dT 



to 



,2r 



(g>(w)) 2 J t0 {l-b{r)f'{w)f 



f'"H + 



b(t) 



i 



2a 2 (^H) 2 J Hto) (l-xf'(w)Y 



2a 2 (^H) 2 A(t ) (i~^/'H) 2 
, 3(/"M) 2 1 



f'H + 



3(r( W )) 2 6(r) 
1 - 6(r)/'(w) 

3(/"H) 2 x 



rfr 



(l-x/'H) 
3(/"H) 2 



dx 



f'H 



f'(w) (1 - xf'(w)) 



dx 



, 3 (f M) 2 



/'"H 3(/"W) ! 



2 (/'H) 2 \(i~xf(w)r 



j» (/'H) 2 
1 

- 1 

3 
~ 2 



(i-x/'H) 



&(*)/"» 



6(to)/'"H 



b(t )f"(w 



b(t ) 

b(t)f'(w) 
1 - h(t)/'H 

2\ 



i-&(*o)/'H 2\i-i(td)/'M. 
(-(5 5t )H + (^ )H) 



2a 2 (^H) 2 
3s" 



where we used the identity (14ot in the last step. Hence 



II 



2a^ 



to 



M zzz (r, z 3 {r)) \dr = ~Y J *j) 



3 = 1 



To calculate integral III, put c 



f( w i)-f( w i) ' 



(47) 

Then for each j ^ I we get a 
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contribution to the sum in III which takes the form 
-—J t0 SAT)Sl{T) { MrJ-z.Mja 

M z (T,Zj(T))+M z (T,Zl(T)) 

(z 3 (r) - zdrW 



dr 



sjsi f m ( 2cf'(w j )f'(w l )-(f'(w j ) + f'(w l )) 



| 2(cf'( Wj ) - l)(cf( Wl ) - 1) 
[x — c) 3 

Now for any s > we can simplify 

Ms) p c /'too/>o~(/'K) + rM) 

\ (x — c) 2 

2(cf(w J )-l)(cf'(w l )-l) 
(x — c) 3 

1 1 



dx. (48) 



dx 



(2cf(w j )f(w l ) - (/'to) + /'to))) , /i( s) _ r 
+ (c 2 /'K-)/'M - c(/'K) + /'to)) + 1) ( (^^e). 

(2c/'to)/'M - (/'(«;,•) + /'to))) (rf>00 2 - c 2 6( S )) 



c 2 (6(s) - c) 2 

„2 r/ 



+ (c 2 /'to)/'to) - c(/'to) + /'to)) + 1) (2cfc(«) - 6(s) 2 ) 

! (-(a: - b(s)) 2 + c 2 (6( s )/'to) - l)(&( S )/'to) - 1)) 



c 2 (6(s) -c) 



2 



9' s (wj)g' s (wi) 1 



= (/to) - /to)) . , , . 

\{9s{Wj) - 9s(wi)Y to - w ') 
= (/to) - f{m)f{Sg s ){w 3l wi), 

so by equations i|48ll and (14511 . 

m = E ^(J^ M Wih)fa,vn) - (s 9t0 )( Wj , wl )) 

= -^J2 s MSh t t °)(z j ,z l ). (49) 

Inserting these integrals into equation lHf))) gives l(23jl . | 

We are now ready to prove the main result, Theorem l2.3l Let s = (s%, . . . , Sk) € 
C fc , z = (zi, . . . , Zk) G Q k and t = (t%, . . . , t*), where < t k < t k -i . . . < ii, be 
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given. We will prove that the characteristic function 



(j)(t, s, z) = E 



exp < i 



of the random vector U — ((Y tl , 7377)7 ■ • • j (Xt k , ~^~)) ls the characteristic func- 
tion of a Gaussian vector. Since we have assumed that Y is an initial measure 
it follows that 4>q(s,z) = cf>(0,s,z) — exp |i Ylj=i s i(^0) r^j - )}- With the con- 
vention tk+i = 0, define the functions <p? , j = 1, . ..,fc, depending on the 



variables s^' = (s^ 



a) 



and the single time vari- 



able t, < t <tj by the following expression: 



<^)( s O) ;Z 0-)) =e 



fe 1 J 1 
exp li Sm ( Y t™ > — > + * E s ™ 



m— j + 1 



m — 1 



0)' 



Israelsson's Theorem Ij3.1|) states precisely that the satisfy equation l|22l) 
with initial conditions tfit^+i ( s ^> = Yt ^( s ^i s j+i> z j+i) f° r J = 
1, ...,k— 1, and 4>^ {s^ k \ z^) — 4> {s^ k \ z^). Thus we may successively 
integrate k times to obtain <j){t, s, z) — <fy^[si,z-\) in terms of the initial condi- 
tions, using Lemma (13. 2J1 in each step. More explicitly, for j = 1, . . . , k — 1 we 
have by Lemma 13.2|l 



x exp < z 



E 

Z=l 



7 



-5 EE -MAffi}, (50) 

1=1 m=l J 



where — I f 2 _ i\ %' / z '"' ; a „j A(J) _ aO') 

wnere — 2 ^ ^ ana ly lm ~ yi 



13+1 (*«>) 



mi [3 



Applying formula lj5T))) fe — 1 times, starting with s, z) = ^^(si, zi), we 
obtain 



*(*,,,,) = *«) n«p |<e-w -Jee «a21. 



i=i 



i=i 



Z=l m=l 



(?) (?) 

where = Zj, Sj 



(51) 



Sj and 



,«) - I,' 



(^m) — ^tljC^m) 



s ™ nf=m ^t^ 1 ( z ™ ) = s "Xl ( z ™) for to = 1, . . . , j 
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After a final application of Lemma l-'t.2t , equation 15111 becomes 

k / j -, j j 



s, z) = m><(*), Hz)) ex P E y E *W - \ E E 

I j=l \ i=l (=1 m=l 

too^w, /WW) exp|*EE -FVF } - \ E E E -M^w 



2 

i=l j"=Z i = l m=l j=max(l,m) 



(52) 



where s • ftj(z) = (si/i' tl (^i), . . . s fe /4 fe (z fe )) and ft t (z) = {h tl {zi), . . . ,h tk {z k )). 
Now 



1 (1 



\ & for' v%w)) 



KI- 1 ) s 'i( log ' i!; "' ( "' ;(fl,) )' 



(/ io g n^ +i '(^^))j =^iog((^ +i °---°^ +i )'(^))=^iog(^(^)), 



3=1 

we have found that 

k 



E^VF } = ^(|-l) S ^logK(^)). (53) 

To evaluate s/s m A; TO = Y^j= m axQ,,m) s f' s m Ajf we distinguish between two 
cases. First, suppose that ht^Zi) ^ h tm (z m ). In this case, if we assume Z > m, 
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we can write 

k 

sis m Ki m = s/ s^Aj^ 

j—max(l,m) 

k „ 

= Sl Y,h%' {zi)s m h 1 ^ {z m )~ (Sh% +1 ) (h^( Zl ),ht(z m )) 

3=1 

2 / / 

X 9(^(z ; ))9(^ m (z m )) ° S V '" / 

2^ g 2 . ^Ir^f^))-^ 1 ^^)) ^ 
2 a 2 , A/^>H^\ 

2 a 2 . / /jr^wjr^ ^ 

2 / 

= sis m -hl l m (z m )(Sh tl )(z h h t t l m (z m )), 

as claimed. Secondly, consider the case h tl { z i) — ht m {z m ). Using the identity 
{S(f o g)) (z) = {g'(z)f (Sf) (g(z)) + (Sg) (z) for the Schwarzian derivative of 
a composition, we have 

k 

SlS m Al m = S l^ S rnMm 
j—max(l,m) 

3=1 

3=1 
k 

= ^lE h t'(^)l ({shir) w (sh%) ( Zl )) 

3=1 

2 * ' 

= sis m -h t ' m (z m ) (Sh tl ) (zi,zi). 
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Thus equation (lo2fl can be expressed 
<f>(t,s,z) 



exp | 1 E s i ) ( Y ° ' . _ h f 

k k | 

^EE^f^-'^) (^t,) (^fe)) | , (54) 



which shows that U is Gaussian with mean and covariance as claimed. 
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